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U INTRODUCTION 


The computer program PLDRVR implements the Hodges-Artley Spring Small Grains 
Planting Date Distribution Model (ref. 1). This docuitient supplements the 
user’s guide (ref. 2) with detailed information with respect tu the PLDRVR 
program. Sample input and output files are available in reference 2. 
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2. PROGRAM DESCRIPTION 


The PLDRVR program is written in the Fortran programni ng language. It was 
developed at the NASA Johnson Space Center on the Earth Observations Division 
Laboratory Systetn AS/3000 computer. The system has an IBM Fortran-H compiler 
and the Conversational Monitoring System (CMS). 

Figure 1 illustrates the arrangement of the program components. The program 
has four iiuijor parts: input (METRDR and SPIN); parameter calculation (SPONGE, 

DEGDAY, and SPFILL); planting nwdel (PLANTM and its subprogram family); and 
output (ANSWER and, optionally, TEST). The iiain program calls each of these 
subroutines once during each iteration. The program stops when the end of the 
ineteprological data file is reached. 

Three intrinsic functions are utilized: MOO, the integer remaindering func- 
tion, in SPLIT and MOSES; EXP, the exponential function, in VAPOR; and FLOAT, 
the integer- to-floating point conversion function, in MOSES. 

Input is from two disk files which are defined in CMS as units 1^ and 20. 

These files contain initial values for the sponge variable and meteorological 
data, respectively. Output is to two disk files which are defined in CMS as 
units 7 and 8. File 7 contains a year's list of the meteorological variables 
and model paraineters to aid the verification of the implementation. 

The PLDRVR program requires a virtual machine with less than one megabyte of 
storage. The amount of temporary disk space required to execute the program 
depends mainly on the size of the meteorological data set. 
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Figure 1.- Hierarchy of Subprograias In the PLDRVR prograa. Dashed/box denotes an optional 
subroutine. Huid)ers refer to the sections In which each subprogram Is discussed. 



3. COMPLETE SUBPROGRAM INFORMATION 


A complete description of each subprogram In the PLDRVR program Is Included in 
this section. It iihould be noted that all variables are single precision and 
that there are no coiniwn blocks. 


3J. MAIN 


The purpose of tha MAIN program Is to access the Input, parameter calculating, 
model, and output subroutines for a location, and to return to the top of the 
program and repeat the process for each location in the ineteorological data 
file. 


CALLING PROCEDURE; Not applicable. 

INPUT PARAMETERS; Not applicable. 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BY; Not applicable. 

SUBPROGRAMS REFERENCED ; Subrrutines ANSWER, DEGDAY, METRDR, PLANTM, SPFILL, 
SPIN, SPONGE, and, optionally, TEST. 

INPUT/OUTPUT DEVICES: None. 

ARRAYS; MDL(IOO), PARMS(4,366) , PLANT(22), and WX(6,366). 

LOCAL VARIABLES ; The following abbreviations are used in tables throughout 
this document: 

A s* alphanumeric 
I * integer 
R * real 

Local variables for the MAIN program are as follows; 

Na>ne Type Description 

CAP R Water-holding capacity of the sponge (moisture variable). 
Initialized in a data statement. 

CRD I Crop reporting district of the segment. 
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Name 

DIV 

IPO 

ISEG 

LAT 

LPO 

MOL 

MPD 

NUMPD 

PARMS 


PINDEX 

PLANT 

SPINIT 

ST 

WX 


Type Description 

A State cliiDiilological division of the weather station. Declared 
as an integer. 

I Initial planting date found by the model. 

I Segment identification number. 

R Latitude of the weather station. 

I Last planting date found by the model. 

I 100-element array of the planting days found by the model. All 

elements initialized to 0 in a data statement. 

I Median planting date found by the model. 

I The total number of planting dates found by the model. 

R 4-by-366 array of parameters used by the planting iiKlel. For 

PARMS (I,J), 0 is the day of year and I refers to a particular 
variable; 1, sponge-precipitation variable; 2, drought 
threshold for range- adjusted base 32®F growing degree days; 3» 
base 32®F growing degree days; and 4, range-adjusted base 32®F 
growing degree days. All elements are initialized to 9999 in a 
data statement. 

I Number of representative dates found by the model. 

I 22-element array of the representative dates found by the model. 

R Starting value for the sponge. 

I State identification number for the segment. 

R 6-by-366 array of daily meteorological variables. For WX(I,J), 

J is the day of the year and I refers to a particular weather 
variable: 1, maximum temperature, ®F; 2, minimum temperature, 

®F; 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. All 
elements initialized to 9999 in a data statement. 

I Year of the weather data. 
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3.2 METRDR 


The METRDR subroutine reads meteorological da*-’\ for one location and year and 
puts them into the WX array, and enu*.; the program after the last set of 
meteorological data has been read and run through the planting model. 

CALLING PROCEDURE ; CALL METRDR(ISEG.WX,LAT,OIV,YR,ST,CRD) 

INPUT PARAMETERS; None. 

OUTPUT PARAMETERS ; CRD, DIV, ISEG, LAT, ST, WX, and YR 

REFERENCED BY ; The METRDR subroutine is referenced by the MAIN program. 

SUBPROGRAMS REFERENCED; None. 

INPUT/OUTPUT DEVICES : Unit 20. file of meteorological data for input. 

ARRAYS ; P(16), T(2,10), and WX(6,366) 

LOCAL VARIABLES ; 

Name Type Description 

CRD I Crop reporting district of the segment. 

DIV A State climatological division of the weather station. 

ISEG I Segment identificaUon number. 

lYR I Year of the weather data. Read from the data and control cards on 

unit 20. 

J I DO loop index. 

JD 1 Day of year. 

K I Loop index. 

LAT R Latitude of weather station. 
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Natne Typ e 


Description 


NCARDS I 
ND I 

P I 

ST I 

T I 


WX R 


YR I 


Loop counter of the nuitiber of data cards read from unit 20, 
Day-of-year counter. 

l6-e1ement array for temporary holding of precipitation data. 

State identification number for the segment. 

2-by-lO array for temporary holding of maximum and minimum 
temperatures. For T(I,J), J is the day of the year, I » 1 for 
maximum temperatures, and I * 2 for minimum temperatures. 

6-by-366 array of daily meterological variables. For WX (I,J), 
J is the day of the year and I refers to a particular vyeather 
variable: 1, maximum temperature, ®Fi 2, minimum temperature, 

°F', 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 

Year for which weather data was obtained. 
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3.3 SPIN 


The SPIN subroutine initializes the sponge at a given location and year, 
either at half-capacity or at a value read from a file. 

CALLING PROCEDURE ; CALL SPIN(STN,YR,SPINIT,CAP) 

INPUT PARAMETERS ; CAP, STN, YR 

O UTPUT PARAMETERS ; SPINIT 

REFERENCED BY ; The SPIN routine is referenced by the MAIN poqram. 
SUBPROGRAMS REFERENCED : None. 

INPUT/OUTPUT DEVICES : Unit 19, file of locations and values at which to 
initiate the sponge for a specific location. 

ARRAYS : None. 

LOCAL VARIABLES: 


Name 


Description 

CAP 

R 

Water-holding capacity of the sponge. 

ID 

I 

Identification number of the sponge value in file 19. 

SPINIT 

R 

Initial value for the sponge. 

STN 

R 

Location identification number from MAIN program. 

VALUE 

R 

Initial value of the sponge read from file 19. 

YR 

I 

Year for which the weather data were obtained. 


r 
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3.4 SPONGE 


The SPONGE subroutine calculates daily values of the sponge and puts thetn into 
the WX array. 

CALLING PROCEDURE; CALL SPONGE(WX,CAP,SPINIT) 

INPUT PARAMETERS ; CAP, SPINIT. and WX 
OUTPUT PARAMETERS; WX (row 6 only) 

REFERENCED BY ; The SPONGE subroutine is referenced by the MAIN program. 
SUBPROGRAMS REFERENCED; Function EP. 

INPUT/OUTPUT DEVICES; None. 


ARRAYS: 

WX(6,366) 

LOCAL VARIABLES: 

Name 

Type 

Description 

CAP 

R 

Water- hoi ding capacity of the sponge. 

JD 

I 

DO loop index (day of year). 

PPN 

R 

Daily precipitation value, inches. 

SPINIT 

R 

Starting value for the sponge. 

TN 

R 

Daily minimum temperature, ®F. 

TX 

R 

Daily maximum temperature, ®F. 

WX 

R 

6-by-366 array of daily meteorological variables. For WX(I,J), 
J is the day of the year and I refers to a particular weather 
variable: 1, maximum temperature, ®F; 2, minimum temperature, 
®F; 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 
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3.4.1 EP 


The EP function estimates monthly pan evaporation from maximutn and minimum 
temperatures using the Trenchard algorithm. 

C ALLING PROCEDURE : EP(TX,TN) 

INPUT PARAMETERS ; TN, TX 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BY ; The EP function Is referenced by subroutine SPONGE. 

SUBPROGRAMS REFERENCED ; Function VAPOR. 

INPUT/OUTPUT DEVICES ; None. 

ARRAYS; None. 

LOCAL VARIABLES ; 

Name Type Description 

TN R Minimum temperature in ®F. 

TX R Maximum temperature In ®F. 
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3.4.2 MPOR 


The purpose of the VAPOR function is to calculate the saturation vapor 
pressure over water in millibars at a temperature in degrees Fahrenheit. 

CALLING PROCEDURE ; VAPOR (T) 

INPUT PARAMETERS; T 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BY : The VAPOR function is referenced by function EP. 

SUBPROGRAMS REFERENCED ; Machine function EXP. 

INPUT/QUTPUT DEVICES : None. 

ARRAYS; None. 

LOCAL VARIABLES ; 

Naiae Type Description 
T R Temperature in ®F. 
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3.5 DEGDAY 


The purpose of the DEGDAY subroutine Is to accumulate growing degree days and 
range-adjusted growing degree days from a starting date with a specified base 
temperature and put the daily values into the PARMS array. 

CALLING PROCEDURE ; CALL DEGDAY (WX, START, BASE, PARMS) 

INPUT PARAMETERS ; BASE, START, AND WX 

OUTPUT PARAMETERS ; PARMS (rows 3 and 4) 

REFERENCED BY ; The DEGDAY subprogram is referenced by the MAIN program. 
SUBPROGRAMS REFERENCED : None. 

INPUT/OUTPUT DEVICES; None. 

ARRAYS; PARMS(4,366) and WX(6,366) 

LOCAL VARIABLES ; 

Name Type Description 

BASE R Base temperature of the growing degree days. 

DD R One day's degree day value. 

I I DO loop index (day of year). 

PARMS R 4-by-366 array of parameters used by the planting nwdel , For 

PARMS(I,J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought threshold 

for range- adjusted base 32®F growing degree days; 3, base 32®F 
growing degree days; and 4, range- adjusted base 32®F growing 
degree days. 

R2GDD R One day's range-adjusted degree day value. 
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Name Type 


Description 


START I 
TM R 

TN R 

TX R 

WX R 


Y R 

I R 


Day of year on which to begin degree day accumulation. 

Daily mean temperature. 

Daily minimum temperature. 

Daily nwximum temperature* 

6-by-366 array of daily meteorological variables. For WX(I,J), J 
is the day of the year and I refers to a particular weather vari- 
able: 1, nwximum temperature, °F; 2, minimum temperature, "Fi 
3, precipitation (liquid equivalent), inches; 4, defined by user; 
5, defined by user; and 6, sponge value, inches. 

Cumulator of growing degree days. 

Cumulator of range- adjusted growing degree days. 
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3.6 SPFILL 


The purpose of subroutine SPFILL is to calculate the drought threshold for the 
range-adjusted growing degree days and the sporfge-precipitation variable and 
put their daily values into the PARMS array. 

CALLING PROCEDURE ; CALL SPFILL (WX,P\RMS) 

INPUT PARAMETERS ; WX 

OUTPUT PARAMETERS : PARMS (rows 1 and 2) 

REFERENCED BY ; The SPFILL subprogram is referenced by the MAIN program. 
SUBPROGRAMS REFERENCED : None. 

INPUT/OUTPUT DEVICES ; None. 

ARRAYS; PARMS(4,366) and WX(6,366) 

LOCAL VARIABLES ; 

Name Type Description 

DELT I Number of days since last increase or saturation in the 
sponge. Current day counts as one day. 

I I DO loop index (day of year). 

IB I First date of sponge wetting or saturation period. 

IDAY I Day-of-year counter in determination or LASTSP. 

INDEX I DO loop index (day of year). 

LASTSP I Last day of sponge wetting or saturation period. 
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Name Type Description 

PARMS R 4-by-366 array of parameters used by the planting nwdel. For 

FARMS ( I, J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought 
threshold for range- adjusted base 32 "F growing degree days; 

3, base 32®F growing degree days; and 4, range- adjusted base 
32“F growing degree days. 

PRE R Sutn of daily precipitation during the period of in which the 

sponge is increased or saturated. 

WX R 6-by-366 array of daily meteorological variables. For WX(I,«1), 

J is the day of the year and I refers to a particular weather 
variable; 1, maximum temperature, °F; 2, minimun temperature, 
®F; 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 
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3.7 PLANTM 


The purpose of the PLANTM subroutine Is to serve as a driver for the planting 
model subprograms. 

CALLING PROCEDURE ; CALL PLANTM ( WX , P ARMS , I PD ,MPD . LPD , PLANT . P I NDEX ,M0L , NUMPD ) 
INPUT PARAMETERS ; PARMS and WX 

OUTPUT PARAMETERS ; IPD, LPD, MOL, MPD, NUMPD, PINDEX, and PLANT 

REFERENCED BY ; The PLANTM subroutine is referenced by the MAIN program. 

SUBPROGRAMS REFERENCED ; Integer functions FSTPLT, LSTPLT, and MODLST and 
subroutines MODMED and PLANTR. 

INPUT/OUTPUT DEVICES ; None. 

ARRAYS ; MDL(IOO), PARMS(4,366) , PLANT(22), and WX(6,36f) 


LOCAL VARIABLES: 


Name 

IPD 

LPD 

MDL 

MPD 

N 


Type Description 

I Initial planting date found by the model (original and final 
estimates). 

I Last planting date found by the model. 

I lOO-element array of the planting days found by the model. 

I Median planting date found by the model. 

I Number of modeled planting days between the first planting date 
and the median date. 


NEWPD I Revised estimate of the initial planting date (includes the 
check for drought). 
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Naii^e Type Description 

NUMPD I Initially, the nunter of planting udtes between the median and 

the last planting dates (12 in the nondrought situation). It is 
then added with N to give the total number of planting dates 
found by the model . 

FARMS R 4-by-366 array of parameters used by the planting model. For 

PARMS(I,J), 0 is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought thresh- 
old for range- adjusted base 32®F growing degree days; 3, base 
32®F growing degree days; and 4, range-adjusted base 32®F 
growing degree days. 

PINDEX I Number of representative dates found by the model. 

PLANT I 22-element array of the representative dates found by the model. 

WX R 6-by-366 array of daily n^eteorological variables. For WX(I,J), 

J is the day of t\ year, and I refers to a particular weather 
variable: 1, nwximum temperature, ®F; 2, minimum temperature, 

“F; 3, precipitation (liquid equivalent). Inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 
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3.7.1 FSTPLT 


The purpose of the integer function FSTPLT is to find the first planting date 
by either of two algorithms (with or without the drought adaptation). 

CALLING PROCEDURE ! FSTPLT ( WX, PARMS , I TYPE ) 

INPUT PARAMETERS ; ITYPE, PARMS, and WX 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BV> The FSTPLT integer function is referenced by the PLANTM 
subprogram. 

SUBPROGRAMS REFERENCED : None. 

INPUT/OUTPUT DEVICES : None. 

ARRAYS: PARMS (4,366) and WX (6,366) 

LOCAL VARIABLES : 

Name Type Description 

ITYPE I Flag for algorithm type: 0, original algorithm (ref. l)j 1, 

modified algorithm for drought. 

JD I DO loop index (day of year). 

PARMS R 4-by-366 array of parameters used by the planting model. For 

PARMS(I,J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought threshold 

for range-adjusted base 32®F growing degree days; 3, base 32®F 
growing degree days; and 4, range- adjusted base 32“F growing 
degree days. 
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Name 


Type Descri ption 

WX R 6-by-366 array of dally meteorological variables. For WX(I,J), 0 
is the day of the year and I refers to a particular weather 
variable: 1, maximum temperature* 2, minimum temperature, 

‘‘Fi 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 0, sponge value, inches. 


3-17 


3.7.2 


MODMEO 


The purpose of the MODMED subroutine is to determine the number of planting 
days between the initial and median planting days and thus obtain the median 
planting date. 

CALLING PROCEDURE ; CALL MODMED (WX, FARMS, I PD. NEWPD»MPO,N) 

INPUT PARAMETERS ; IPD, NEWPD, FARMS, and WX 
JUTPUT PARAMETERS ; MPD and N 

REFERENCED BY ; The MODMED subroutine is referenced by subroutine PLANTM. 
SUBPROGRAMS REFERENCED; Function MEOPLT. 

INPUT/OUTPUT DEVICES ; None. 

ARRAYS ; PARMS(4.366) and WX(6,366) 

LOCAL VARIABLES; 


Name 

Jm. 

hescription 

IPD 

I 

Initial planting date found by the model without a check for 
drought. 

MPD 

I 

Median planting date found by the model. 

N 

I 

Time in modeled planting dates between initial and median 
planting dates. 

NEWPD 

I 

Revised estimate of the initial planting date (includes check for 
drought). 
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Name Type Description 

PARMS R 4-by-366 array of parafneters used by the planting model* For 
PARMS(I,J), J Is the day of year and I refers to a particular 
variable: 1, sponge-participation variable; 2 , drought threshold 
for range- adjusted base 32®F growing degree days; 3, base 32®F 
growing degree days; and 4, range- adjusted base 32®F growing 
degree days. 

WX R 6-by-366 array of dally meteorological variables. For WX(I,J), J 
Is the day of the year and ! refers to a particular weather 
variable: 1, maximum temperature, ®F; 2, minimum temperature, 

®F; 3, precipitation (liquid equivalent). Inches; 4, defined by 
user; 5, defined by user; and 6, sponge value. Inches. 
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3. 7. 2.1 MEDPLT 


The purpose of the MEDPLT function is to established the median planting date 

as occurring a certain number of days after the first planting date. 

CALLING PROCEDURE ; MEOPLT ( WX . PARMS . I PD , N ) 

INPUT PARAMETERS ; I PD. N, PARMS. and WX 

OUTPUT PARAMETERS; Not applicable. 

REFERENCED BY ; The MEDPLT function is referenced by subroutine MODMED. 

SUBPROGRAMS REFERENCED; None. 

INPUT/OUTPUT DEVICES ; None. 

ARRAYS; PARMS(4,366) and WX(6.366) 

LOCAL VARIABLES ; 

Name Type Description 

I I Counter of modeled planting dates. 

IPD I Initial planting date found by the iiwdel . 

N I Number of planting days between the initial and median planting 

days. 

PARMS R 4-by-366 array of parameters used by the planting model. For 

PARMS (I.J). J is the day of the year and I refers to a particular 
variable; 1. sponge-precipitation variablej 2, drought threshold 
for range- adjusted base 32®F growing degree days; 3, base 32®F 
growing degree days; and 4, range- adjusted base 32“F growing 
degree days. 
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Name Type Description 

WX R 6-by-366 array of daily meteorological variables. For WX(I,J), J 
is the day of the year and I refers to a particular weather 
variable: 1, maximum temperature, “Fi 2, minimum temperature, 

®F; 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defimjd by user; and 6, sponge value, inches. 


3-21 


3.7.3 LSTPLT 


The purpose of integer function LSTPLT is to establish the last planting date 
as occurring a certain nuinber of days after the median. 

CALLING PROCEDURE; LSTPLT ( IF IRST , WX , PARMS , NUMPD ) 

INPUT PARAMETERS ; IFIRST, NUMPD. PARMS, and WX 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BY ; The LSTPLT integer function is referenced by subroutine 
PLANTM. 


SUBPROGRAMS REFERENCED; None. 
INPUT/OUTPUT DEVICES ; None. 

ARRAYS: PARMS(4,366) and WX(6,366) 


LOCAL VARIABLES : 

Name Type Description 

IFIRST I Day from which to locate the last planting date. 

JD I GO TO loop index (day of year). 

N I Counter of modeled planting dates. 

NUMPD I The number of planting dates between the median and final 

planting dates. 

PARMS R 4-by-366 array of parameters used by the planting model. For 

•MRMS(I,J), J is the day of the year and I refers to a particu- 
lar variable: 1, sponge-precipitation variable; 2, drought 

threshold for range- adjusted base 32®F growing days; 3, base 
32®F growing degree days; and 4, range-adjusted base 32**F 
growing degree days. 
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Descri ption 

6-by-366 array of daily n«teorologica1 variables. For WX(I,J), 
J is the day of the year and I refers to a particular weather 
variable: 1, maximiim temperature, °F; 2, minimum temperature, 

3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 


3.7.4 MODLST 


The purpose of integer function MODLST is to check for drought on the last 
planting day, and, if there is a drought, to recalculate the last planting 
day. 

CALLING PROCEDURE; MODLST ( WX , FARMS ,LDAY ,NUMPD ) 

INPUT PARAMETERS ; LDAY, NUMPD, FARMS, and WX 
OUTPUT PARAMETERS ; NUMPD 

REFERENCED BY; The MODLST integer function is referenced by subroutine 
PLANTM. 


SUBPROGRAMS REFERENCED; None. 


INPUT/OUTPUT DEVICES ; None. 

ARRAYS: PARMS(4,366) and WX(6,366) 


LOCAL VARIABLES: 


Name Type 


Description 


FLAG I 

I I 

JO I 

LDAY I 

NUMPD I 


Precipitation event flag. 

DO loop index (day of year). 

GO TO loop index (day of year). 

Initial estimate of the last planting date. 

The number of planting dates between the median and final dates. 


Name Type Description 

FARMS R 4-by-366 array of parameters used by the planting model. For 
FARMS J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought threshold 

for range-adjusted base 32®F growing degree days; 3, base 32®F 
growing degree days; and 4, range- adjusted base 32°F growing 
degree days. 

WX R 6-by-366 array of daily meteorological variables. For WX(I,J), J 
is the day of the year and I refers to a particular weather vari- 
able: 1, maximum temperature, “F; 2, minimum temperature, ®F; 

3, precipitation (liquid equivalent), inches; 4, defined by user; 
5, defined by user; and 6, sponge value, inches. 
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3.7.5 PLANTR 


The purpose of the PLANTR subroutine is to determine the series (strings) of 
consecutive planting days within the planting period and to obtain dates which 
represent the strings and the planting period. 

CALLING PROCEDURE; CALL PLANTR(WX,PARMS,ISTART, MEDIAN, PLANT, PINDEX,NUMPD,MDL) 
INPUT PARAMETERS ; ISTART, MEDIAN. NUMPD, PARMS, and WX 
OUTPUT PARAMETERS ; MDL, PINOEX, and PLANT 

REFERENCED BY ; The PLANTR subroutine is referenced by subroutine PLANTM. 
SUBPROGRAMS REFERENCED : Subroutines ILAST and SPLIT. 

INPUT/OUTPUT DEVICES : None. 

ARRAYS ; IFIRST(22), LAST(22). LENGTH(22), MDL(IOO). PARMS(4,366) , PLANT(22), 
and WX (6,366) 

LOCAL VARIABLES ; 

Name Type Description 

I I DO loop index and index of modeled planting date strings. 

IFIRST I 22-element array of first dates in each string of planting 

dates. 

IIIPD I Third representative date for a string of planting dates. 

II PD I Second representative date for a string of planting dates. 

INI I Length of time (days) between a solitary planting date and the 

end of the previous string of planting dates. 

IN2 I Length of time (days) between a solitary planting date and the 
beginning of the next string of planting dates. 
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Name Type 


Description 


IPD I 

ISTART I 

LAST I 

LENGTH I 

LTOTAL I 

MDL I 

MEDIAN I 

NSTRNG I 

NUMPO I 

PARMS R 


PINDEX I 
PLANT I 

WX R 


First representative date for a string of planting dates. 

Initial planting date found by the model. 

22-element array of last dates in each string of planting dates. 

22-element array of the number of days in each string of 
planting dates. 

Total number of planting days [the sum of the length (I)'s]. 
100-element array of the planting days fo.'id by the inodel. 

Median planting date found the model. 

Number of strings of planting dates. 

The total number of planting dates to be found by the model. 

4-by-366 array of parameters used by the planting model. For 
PARMS ( I, J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2» drought 

threshold for range- adjusted base 32®F growing degree days; 

3, base 32®F growing degree days; and 4, range-adjusted base 
32®F growing degree days. 

Number of representative dates found by the model. 

22-elen»ent array of the representative dates found by the model. 

6-by-366 array of daily meteorological variables. For WX(I,J), 

J is the day of the year and I refers to a particular weather 
variable: 1, maximum temperature, ®F; 2, minimum temperature, 

°F; 3, precipitation (liquid equivalent), inches; 4, defined by 
user; 5, defined by user; and 6, sponge value, inches. 
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3. 7.5.1 HAST 


The purpose of the HAST subroutine is to determine the last day in a series 
(string) of consecutive planting days. 

CALLING PROCEDURE ; CALL ILAST( IFIRST,LAST,N,WX,PARMS,LTOTAL,NUMPD.MDL) 

INPUT PARAMETERS: IFIRST, LTOTAL, NUMPD, PARMS, and WX 

OUTPUT PARAMETERS ; LAST, MDL, and N 

REFERENCED BY ; The HAST subroutine is referenced by subroutine PLANTR. 
SUBPROGRAMS REFERENCED; None. 

INPUT/OUTPUT DEVICES; None. 

ARRAYS ; MDL (100). PARMS (4, 366 ) , and WX(6,366) 

LOCAL VARIABLES; 


Name 

Type 

Description 

IFIRST 

I 

Day from which to locate the last planting date. 

JD 

I 

GO TO loop index (day of year). 

LAST 

1 

Last date in a string of pi anting dates. 

LTOTAL 

I 

Total number of planting dates. 

MDL 

I 

lOO-element array of the planting days found by the model. 

N 

I 

Length of planting date string in days. 

NUMPD 

I 

The total number of planting dates to be found by the model. 
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Description 

4-by-366 array of parameters used by the planting model. For 
FARMS d is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought threshold 
for range- adjusted base 32 ®F growing degree days; 3, base 32 “F 
growing degree days; and 4, range- adjusted base 3?®F growing 
degree days. 

6-by-366 array of daily meteorological variables. For WX(I,J), 0 
is the day of the year and I refers to a particular weather vari- 
able; 1, maximum temperature, ®F; 2, minimum temperature, ®F; 

3, precipitation (liquid equivalent), inches; 4, defined by user; 
5, defined by user; and 6, sponge value, inches. 



3. 7. 5. 2 SPLIT 


The purpose of subroutine SPLIT is to split a series (string) of consecutive 
planting days into substrings based on the overall length of the string* and 
to obtain dates which represent these substrings of days. 

CALLING PROCEDURE; CALL SPLIT(IF.MED,IL,N,DAY1,DAY2,DAY3) 

INPUT PARAMETERS ; IF* IL* MED, and N 

OUTPUT PARAMETERS ; DAYl, DAY2. and DAY3 

REFERENCED BY ; The SPLIT subroutine is referenced by subroutine PLANTR. 
SUBPROGRAMS REFERENCED; Function MOSES and inachine function MOD. 
INPUT/OUTPUT DEVICES; None. 

ARRAYS; None. 

LOCAL VARIABLES; 

Name Type Description 

DAYl I First representative date for a string of planting dates. 

DAY2 I Second representative date for a string of planting dates. 

DAYS I Third representative date for a string of planting dates. 

IF I First day in first substring of planting days. 

IIF I First day in second substring of planting days. 

II IF I First day in third substring of planting days. 

IL I Last day in a string of planting days. 

K I Length of first substring of planting days. 

KK I Length of second substring of planting days. 
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Name Type Description 

KKK I Length of third substring of planting days. 

KKMOD I Flag for odd or even length of substrings: 0 if even, nonzero if 

odd. 

MED I Overall median planting date found by the model. 

N I Length of planting date string in days. 
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3. 7. 5.3 MOSES 


The purpose of the MOSES function is to determine the median of a series 
(string) of consecutive days. In the event of a string of even length, the 
median is rounded toward the overall median. 

CALLING PROCEDURE; MOSES (N . I I .MED , IL ) 

INPUT PARAMETERS ; II. IL. MED, and N 

OUTPUT PARAMETERS ; Not applicable. 

REFERENCED BY; The MOSES function is referenced by subroutine SPLIT. 
SUBPROGRAMS REFERENCED: Machine functions MOD and FLOAT. 

INPUT/OUTPUT DEVICES: None. 


ARRAYS; None. 
LOCAL VARIABLES; 


Name 

Im. 

Case riot ion 

II 

I 

First day in substring. 

IL 

I 

Last day in substring. 

IREMDR 

I 

Flag for odd or even length of substring; 0 if even, nonzero if 
odd. 

MED 

I 

The overall median planting date found by the model. 

N 

I 

Length of planting date string in days. 
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3.8 ANSWER 


The purpose of subroutine ANSWER is to write the results of the niodel for one 
location to a disk file. 

CALLING PROCEDURE ; CALL ANSWER(YR.ISEG,IPD.MPD,LPD, PLANT, PINOEX, MOL. NUMPD) 
INPUT PARAMETERS ; IPD, SEG, LPD, MOL. MPO, NUMPO, PINDEX. PLANT, and YR 
OUTPUT PARAt-iETERS ; None. 

REFERENCED BY ; The ANSWER subroutine is referenced by the MAIN program. 
SUBPROGRAMS REFERENCED; None. 

INPUT/OUTPUT DEVICES: Unit 7. 

ARRAYS : MDL(IOO). PARMS(4,366) , and PLANT(22) 

LOCAL VARIABLES: 


Name 

Type 

Description 

I 

I 

DO loop index (day of year). 

IPO 

I 

Initial planting date found by the model. 

I SEG 

I 

Segment identiflcatio»^ number. 

LPD 

I 

Last planting date found by the model. 

MDL 

I 

100-element array of the planting days found by the model. 

MPD 

I 

Median planting date found by the model. 

NUMPD 

I 

The total number of planting dates found by the model. 
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Name Type Description 

PARMS R 4-by-366 array of parameters used by the planting model. For 

PARMS(1,J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought 

threshold for range- adjusted base 32 “F growing degree days; 3, 
base 32®F growing degree days; and 4, range-adjusted base 32®F 
growing degree days. 

PINDEX I Number of representative dates found by the model. 

PLANT I 22-element array of the representative dates found by the model. 

YR I Year for which the weather data were obtained. 
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3.9 TEST 


The purpose of the TEST subroutine is to aid the checking of program 
implementation by writing the daily values of WX and PARMS arrays to a disk 
file. 

CALLING PROCEDURE ; CALL TEST (WX, PARMS. I SEG) 

INPUT PARAMETERS ; ISEG, PARMS, and WX 
OUTPUT PARAMETERS : None. 

REFERENCED BY ; The TEST subroutine is referenced by the MAIN program. 
SUBPROGRAMS REFERENCED: None. 

INPUT/OUTPUT DEVICES : Unit 8. 

ARRAYS ; PARMS(4,366) and WX(6,366) 

LOCAL VARIABLES : 

Name Type Description 

I I DO loop index (day of year). 

ISEG I Segment identification number. 

J I DO loop index. 

PARMS R 4-by-366 array of parameters used by the planting model. For 

PARMS(I,J), J is the day of year and I refers to a particular 
variable: 1, sponge-precipitation variable; 2, drought threshold 

for range- adjusted base 32 “F growing degree days; 3, base 32 ®F 
growing degree days; and 4, range-adjusted base 32*’F growing 
degree days. 
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Name Type Description 

WX R 6-by-366 array of daily meteorological variables. For WX(I,J). J 
is the day of the year and I refers to a particular weather vari- 
able: 1, maximum temperature, ®F; 2, minimum temperature, ®F; 

3, precipitation (liquid equivalent), inches; 4, defined by user; 
5, defined by user; and 6, sponge value, inches. 
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APPENDIX 


PLDRVR FORTRAN CODE 



o r» r* -I rtrt n.n on>r>r>r»r» rwi nno n n r»n oaoormriortoo or»r»of»f>nor»of>nn 


FILFI Pl.OffVP FOPTpan 


ORlGlNftL PAGE B 
OF POOR QUALIT» 


C0MVF.PSATI0N4L MONITOR STSTfiM 


OQ069AM impipmfnts thf mooofs-artlct algorithm for spring small 

GRAINS PLANlINfi niSTRIRHTION 

RPAL W< (G» 16^) *L AT *P4RMSI4*.T6A,j 

WX(I*'jV WHPPP J IS JULIAN. 0 ATE ANO I IS AS _ FOLLOWS: 

laT, mAIIa*(IM TE^^Pt'P ATURF 

|«K mTNIviim t|mpe»ATURF ••wx is FIlLFU'IN suhhoutinks 

laS* PPeCIPlTAf ION ••METRDR ANO SPONGE __ _ 

la4 « OoFn •••♦••••••*••••••••••••#•••*•< 

|>s« OPFN 

farms n *J?^S5FHfp‘ j TS JttLlAM OATF ANO I IS AS FOLLOWS: 

1*Tf SPONGP/PPECtPlTATinN VARIABLE, , 
l«2. OOOtjGHT Gl)')32R THRESMOLO. 45*SPQN0E 

ORppTNG.nEOeEt’S OAYS. <5^^^ _ 


:aM OATF ANO 
‘ATION VaRIA 


AS follows: 


mu, PAN(iF-AnJU*;TEO GOOWINGE OFGREE OAVS (G0032R) ♦ BASE 32 
IS FHLFO alTMSlMROUTINF.S f3FGOAy ANPsPFILL 


PARMS is FULFO alTMSlMROUTINF.S f3FGOAy ANPSPFILL 
JNTFGFP yP*PI.ANT (22) .PIMnEX^STtCROfMOL ( 100) 

DATA HOL/)nO*0/.PARMS/l464*R9R9,0/ 

OATA WX/2)9A*9999.0/ 

CAp‘’lS*THF.^TofAL MOISTUPE CAPACITY OF THE SPONGE 


INPUT 


19 SPONGE values— INIT IAL^VALUES OF SPONGE, USUALLY 
THE VALUE ON THE ENO OF. THE, PREVIOUS YEAR, OTHERWISE 
THE SPONGE IS initialized AT HALF-FULL. VALUES ARE 
PY SEGMENT. _ 


OUTP)iT file* 


FILF 20 year met — YFaR»S WEATHER DATA BY SEGMENT 
FTCE 7 --THF OUTPUT FROM HODGES-ARTLEY SPRING 
GRAINS PLANTING ALGORITHM 


' ?ALL*MFfROO(rSEG.WX*LAT.niV.YR.STfCRO) 

FILL PART OF WX. ALL VAR I ARLES ARE OUTPUT^^ 
T8FGxSEGMF^»T NIIMREP. LAT*STATfnN LATITUDE. YR«YEAP 
•CALL SPIN( ISEG.yP.SPINiT.CAP) 

GET initial value EOR ThE SPONGE (SPINIT) FOR THIS YEAR 


CALL SPONGE (WX.CAP.SPINIT) 

PUT daily va( (ifs of sponge into WX 

- CALL 0EG0AY(WY.lO2.f),PA0MS) , . . , ^ 

CALCULATE CUMULATIVE GROWING DFGREE DAYS FROM DAY 1. WITH A BASE 
OF 3? OEG. F. AMO PUT ANSWERS IN PARMS, 00 FOR REGULAR ANO RANGE 
ADJUSTED GoOwimG OFGPEE DAYS. 

■“CALL SPEIl L (V^ .PARMS) 

Fill PAPMS ifTM THE daily SPONGE/PRECIPITAT ION VARIABLE ANO DROUGHT 
TMPFSHHOLD OF GDD32R. _ . ... 


CALL PL ANTM(W<.PAHMS. IPO.MPO.LPO. PLANT. P INDEX. MOL. NUMPO), 

PUN TMF PLANTING MODEL WITH WX AND PARMS ELEMENTS. OUTPUT IS 
IPD- INITIAL PLANTING DATE. mpq-meDIAN PLANTING DATE. LPO-LAST 
planting OATF. PLANT-OATES WHICH REPRESENT THE MODELEO PERIOD. 
PINOFx-NUMPFo OE RFPRESFNTaTIVE OATES. MOL-EaCH MODELED PLANTING DAY 
IN TWE PERIOD. MUMPO-NUmbEP OF DAYS IN TH£ PERIOD. 

PASS THESE Values (with ioentieication) to answer for output. 

, call ANSWFP(YR.TSeG.IPO.MPO.LPf). plant. PINo£X.mOL*NUMPO> 
MFPF.SupooUTTrlE TEST MAY BE CALLED SO THAT VALUES IN 

M* ano,paoms can be checked. 

GO TO I 
ENO 


L RXRMS ARRAY*********************************** 

SUBPOUTTNF SPFILL (WX. PAPMS) 

'U*J0^UTINC CALCULATES THE SPONGE/boECIPITaT ION PARAMETER EOR THE 
ACCOPOING TO THE HOOOES-AR TLEY PLANTING DATE ADJUSTMENT ALGORITH 
"" * ■ “ DROUGHT G0032R THRESHMOLO VALUE 


ALSO CALCULATES 


EOR THE YEAR 
algorithm 
(45*SP0NGE) 


STANOAPO MET. AORAY. PARMS IS OUTPUT 
AL,w* (G. !S6) .PAPMS(4.3G6) 

CALCULATE THF DROUGHT G0032R THRESHHOLO VALUES FOR EACH DAY 
no s TNDE* * I. 

,®AR“S (^* TnOEx) * 45.*WX (6. INDEX) 

It^*’.DATLY VALUE OF THE SPONGE/PRECIPITATION VARIABLE 
DO 100 INOFX * 1, 366 

find LASTSP. THE last OAY OF SPONGE INCREASE 

lOAV a I.iOEX . I 
10 IDAY * IDAY - I 

IF (IDAY .LT. I) GO TO IS 
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CUU OCU c; UC U U CUOUU UVU OCU O C’O 


Ptn9VP P09To*M « 




COMVCRSATIONAL MONITOR SYSTEM 




TF (w*(<s.in^Y-n,..LT 

IF (WX<f».IOAT) ,fcQ, 

r.O TO 10 

TOAV * I 

CONTINUE 

LASTS** ■ IITAY 


, WXIfMinAY)) GO TO ?.0 
A.O .ANO. mIX(6*IDAY-1) 


.EO. B.OI GO TO 20 


FIND 10* The FI9ST DATE OF SPONGE GETTING PE9100 
TB * LASTSP 

35 IF <WX(6.1BI- .LT. WX(6«IR-in GO TO 40 
TP « IP - 1 

IF (19 .GT. 0 ) GO TO 35 
IB * I 

40 CONTINUE 

IB * 19 ♦ I 

SIN haily ppfcioitation during the period of Increasing or saturated 

SPONGE, » 

PPE * n. 

do go i*tb.lastsp 

PRF 3 PHF ♦ *»X( 3 fI) 

60 CONTINUF 

OELT^^k^THF^fENGTM Kf*'t?ME*s1ncE THE LAST SPONGE INCREASE (DEFINED 
WITH THF LAST DAY AS DAY 1). 

FARMS ( 1 ♦INDEX ) * WX( 6 fLASTSP)*PRE/nELT^ 

STORF SPVAW \/AlUE FOR TODAY IN PARMS ( 1 flNDEX) 

100 CONTINUE 
RETURN 

END t 

SU 890 UTINF nEGDAY(WXfSTART# 9 ASEtPARMS> 

SUPROUT INF .Cfl^CULATES_GRpWlNG_nEGREES_nAYS 


(GDO) 
DAYS are 


growing degrfe days (Gdor). the OFGPEF 

starting OATF START WITH A PASF TEMPF^ATURE OK 
THF «ET. OflTA IN standard POSITIONS. 

THE DAILY VALUES OF G ')0 ANO GDOP ARE PLACED 
REAL WX ( 6 , 366 ) ♦PARHS( 4 , 366 ) 

integfr start 

YxO, 

Z* 0 . 

00 5 1 x 1,166 

TF(WX ( I , 1 ) .GF.R 99 R. 0 )RETURN 

if(I,lt.start)go to 5 

TXxWXd.T) 

TN*WXC?,I) 

TMx(TX*TN) /?. 

TM-MEAN Temperature, tx-maximum temperature, 
tooay»s degree days are do, always positive 

OOxTM-HASE 

lF(nr),LT.n.)DOxo. 

Y X y ♦ no 

"ARMS ( 3 , 1 ) X Y 

T 00 AY*S degree days WITH THE RANGE ADJUSTMENT 
ARE RPGDD, always POSITIVE 
o?GDDxnD-.Ol*(TX-TN)**a 
IE(P?GDO.LT.O.)R 2 GDDxO 
7 x 7 ,r?CiDO 
PAPMS( 6 ,I)xZ 

s continue 

PFTURN 

ENO 


base, wx must Have 


IN array farms 


tn-minimum temperature. 


( 1 % OF the daily Range) 




anting sU0ROUTINES»****#«***»««*«*««**«****«*«***«*»*« 

9LANTM(WX.PAPMS.IPD,MPD,LP0, PLANT, P INDEX, MOL, NUMPO) 

SMALL ORAINS 

.pl.nt< 22 . 

IPO X FSToLT |WX,R4>»MS,0) 

NON-OROUGMT cituat ton— initial planting date 

NEWPO X FSTPLTIWX.Oi^RMS.)) 

RECALCULATF the iNlTlnL PLANTING OATF— CHECK FOR DROUGHT 

Call modmfduvx^pawhs, iPn,NF'*PD,MPo,N) 
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nr> or» ooo nor> o r>or» oo ooo 


ORIGINAL PAGE IS 

OF POOR QUALITY 


FILFS PLORVP rnuTOAM 4 


c 

g 

c 


CONVERSATIONAL MONITOR SYSTEM 

» 

N IS NUMREP OF days from NEWPO TO 


C 

C 

C 

c 


CALCHLATF mEOIAM PLANTING 0*TE MPQ* 

Mpn. 

NU'^PO T«1 ThE MiiMBEo OF MOOFLEO PLANTING DAYS BETWEEN TmE MEDIAN AND 
LAST MOOFL^-O planting OATES 

NMMPD * 1? • . . 

LOO ■ L^TplT(mPo*.V**PARmS.NHmPD) . 

LPO « MODI ST(W«*PA6rMS*LPntN»lMP()) 

NUMPO IS NOW JFSrr TO TMfe” TOTAL NUMBER OF MOOELEO PLANTING OATES 
NNMPO « Ntl^pO ♦ N 
fPO « NFWOO * 

Call PLANT«(W<.PAPMS.jPn.MPO.PLAMT*PINOEX.NUMPO.MnL» 
obtain days which REPRESENT THF MODELED PLANTING PERIOD. 

RETURN 

IuRroutin^ PL antr ( wx. p arms. 1ST apt .meoi an*plant.p Index ♦numpo.mol) 

REAL PARMS(4.1G6> 

INTFGFP PINOe*♦IFIPSr(^?» fLA„ , . 

find FIRST OflY. then «0RK FORWaRQ TO 
nstpng*niimrep nF stpings 

LENGTH! n*LEMGTH OF ITH STRING 

ingSTd) ■ ISTART 
LTOTAL * n 
NSTPNG * 0 

DO 200 I X I, 2^ 

200 LENGTH! I) * 0 

.1 ■ .1 

FIND FIRST and LAST DAYS AND NUMRER’ OF PLANTING DAYS IN EACH STRING 


•GfP Pt NOE* ♦ IF IPS T!2?» .LAST! 2?) .LENGTH (23} .PLANT! 22) ♦MDLdOO) 
- strings of planting OATES 


250 


CALL ILASTdFIRSTd) .LAST!! 
LTOTAL a LTOTAL ♦ LENGTH(I) 


I) •LENGTH!!) ♦WX»PARMS.LT0TAL*NUMP0*M0L) 


NSTRNG X NSTRNG ♦ I 

1 IF (LTOTAL .EO. NUMPO) 
X I « 1 

[FIPSTd) » LASTd-l) 
»0 TO 250 


GO TO 300 

* 1 


CHECK FOP 
300 


320 I X 

II 


IF (NSTRNG 
I ■ 0 


SOLITARY PLANTING DAYS 
•EO. 1) GO TO 400 


AND ASSIGN TO NEAREST STRING 


I ♦ 1 
(LENGTH! I) 
(1 .Fo. n 


(I .EO. 


.GT. 1) GO TO 3R0 
GO TO 360 
NSTRNG) GO TO 365 


SOLITARY planting OATE IN 

oetepminf.neapfst string 


between TWO STRINGS 


IF 



(IFTPSTd} 
(IFIHSTd) 
360 


.LT. 

, , _ .GT. 

GO TO 

SOLITARY PLANTING DATE 
360 


- LASTd-l) 

IFiRSTd) 

TO 360 

TO 36S ^ ^ ' 

ING nearest to THE MEDIAN 
MEDIAN) GO TO 360 
MEDIAN) GO TO 365 


IFIRSTCdd 
LENGTH! T ♦) ) 
LENGTH d) X 
GO TO 390 


IFIRST(I) 
LENGTH! dd 


AT BEGINNING OF THE SERIES 

♦ I 


SOLITARY PLANTING DATE AT END 


IFlPSTd) , 

* LENGTM(I-l) 
0 


365 LASTd-J) X 
LENGTH d-n 
LENGTH (I) X 
GO TO 390 

THE NF*T STPIMG 

^ IE (I ,LT. NSTRNG) GO TO 320 
400 CONTINUE 

SPLIT strings to obtain REPRESENTATIVE PLANTING OATES 


3?0 


PTNDEX X 0 

DO Sno T s 1. NSTRNG 
IF (LENGTH!!) ,EO. 0) 


GO TO 500 
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FItE* PLORVP FORT9AM h 


pONVERSATIONAL MOMirOR SVSTEM 




CALL RpLTTa*'i 9 sr(n fMFoiANtLASTd) tLASTn>-iFmsT(n*uiPo.nPo. 


J[IPn» 


■ oIMHEX*! 
PLANf (PIMOF. o » fPO 


[f 

f 

f 

F 


? 


c 

c 


500 CONTrwUE 


Ipn ,nE, 
IPO .HE 
tipn 
ripo 


,NF, 


n> PiNoe* * piMPE’* ♦ 1 
0) PL^NTtPlNOFW) * IIPO 
0) PINDEX « PTMOEX ♦ I 

o» PLANT iPi^Joex) « nipo 


PFTUPN 

END 


F I o S_T 


INTEGFP FHNCTTOM 
rs FUNCTF' 


c THIS function FINDS THE INITIAL 


DATE CALC 
FSTPyT_(WA.PAPMS.l 


TTpb* 

PL ANT INC DATE 


T I 0 N 


planting OATFS MODFL. (FIRST day 0MCE'G00'32R EXC 
PONGE ^ 


s 




Is 

.IN. 


rHEN G 0 D 32 P 
AND SPVAR IS 


C 

c 


c 

c 

c 


c 

c 


V 

8 


I TYPE 
IF 




less than 4, 

PAPMSI?* JO) * 

FSTPLT 

WY (6.1G(S) *PAR'4S<4*1G6) 

. FOP original ALOOPTIHM* 

(ITYPE .EO. n 60 TO 4 


FXE< 

LES« 


;eeds 

i THAN 


FOP THE HOOGES-ARTLET 
EDS 180 OR. IF 
ECTED value OF 
equal TO 2.).) 


COR 

OR 


ITYPE«1 FOR algorithm WITH DROUGHT. 


THE NO-MOISTtIPE (ORIGINAL) ALGORITHM 
DO 3 JO > 1*366 

FSTPLT ■ JO ... 

CHECK PLANTING INITIATION CRITERIA— 

T«0. IS THF O003PR THRESHHOLD CONSTANT, 2.0 IS THE SPONGE/PREGIPITA- 
tlON VAPTAPLF THPESHHOLD.. ... 

IF (PAPMS(4,JU) .GC. lAO.OnO .aNO. PAHMS(ltJO) .LE. 2.00) GO TO 7 

3 CONTINUE • i 

RETURN 

algorithm WHICH INCLUDES MOISTURE/DROUGHT CONDITIONS 

4 CONTINUE 

DO 5 JO « 1 . 366 

„ FSTPlT « JO 

CHECK PLANTING INITIATION CRITERIA 
G0032R >* lao and SPVAR <»p 

“ -- G0032R >« 4S*SP0NGE AND SPVAR <■ 2. 

4,00) .and. 

180,000 .AND, PARMSd.JO) 

- 4,00) .AND. 

& (PARMS(4.J0) ,5E. farms (2, JO) .and. PARMSd.JO) 
h _ 60 TO 7 


OR, IF SPONGF < 4. 
IF ((MX(6,jn) ,r»E. 

(PARMS(4,jr)) .gE. 
IF (IWX(6.J0) ,LT. 
(PARMS(4.J0) ,ge. 


,LE, 2,00)) GO TO 7 


,LE. 2.00)) 


? 


8 

8 

C 

C 


c 

c 


c 

c 

c 


CONTINUE 

RETURN 

END 


C U L A T I 0 N 


••#•••••••••«•••»• median DATE CAL 

SU8R0UT I NF MOOMED ( WX , P ARMS , I po , MENRO . MPD • N) 

THIS SUBROUTINE COMPUTES THE MEDIAN MODELED PLANTING DATE* 

AND COUNTS TwE NUMPER N OF DAYS pETWEEN TmE INITIAL AND MEDIAN P DAYS 
REAL WX (6.166) .PARMS(4. 366) 

N IS THE NUMBER OF MODELED PLANTING DAYS BETWEEN INITIAL AND MEDIAN 
■ “ DAYS 


moisture SITUATIONS 
IPU) GO TO 100 , 
PLANTING DAYS TO THE 


HOOFLED PLANTING 
N K 10 

N IS 10 IM ample 
IF (NFWQD ,GE. 

FINO THE NfjMBFP OF 
UP TO NEXT DAY) 

N « N ♦ O.S * 0,5*(IPD-NFWPO) 

100 MPD * MEOPLT(WX,PARMS,NEWPO,N)r 
RFTURN 
END 

function MFOPLT(WX.PARMS.IPn.N) 

THIS FUNCTTOm FINOS THE MEDIAN PLANTING DATE 

APTLEV planting date- algorithm, (nth day of 
EOUA t TO ? ONCE THE INITIAL PLANTING DATE IS 
REAL WX (6.166) ,PAR!4S(4, 366) 

MEOPLT s IPO 

CHECK PLANTING OflTE CITERlA 
10 IF (PflOMSf I.mfdpLT) .LE. 2.0) 1*1.1 
IF (I .LT. N) MFDPLT * MEOPLT ♦ I 
IF (I ,lT. N) go to 10 
RETURN 


MEDIAN DATE ESTIMATE (ROUND 


ACCORDING TO THE HODGES- 
SPVAR LESS THAN OR 
DETERMINED) 


A -4 



riLF: PLORVR rnRTRAW A 


OF POOR QUALITY 

CONVERSA'f lONAL MONITOR SYSTEM 


END 


C 

C 


INTEGFO F'lMiUTfON 


T OATF CAL.__, 
LSTPLT( IFIPST.-afRARMS^NOMPO) 
CALCULATFS the cast PtAMTiMO 0 
AFTEO THF 1ST mOOELEO PLANT r^•^, 


A T F 


C A 


C U 


A T I 0 N 


HmSoJ * 

OA^E *s’0CCufi»lN0 NUtiRO PLANTING DAYS 
PAY IFIRST. ' ' 

WX ;PARMS(4t3A6T 

IFIRST - I 


CRITERIA 

• LE. 


RETURN 


N ■ N ♦ I 


a 

c 

a 


a 


RfaL 

JO « 

N ■ 0 

100 JO » jn ♦ i 
CHECK PLANTI^'G OATF 
IF (PAOMSn*J0) 

LSTPLT * JO 
IF IN .r»T, NUMPOI 
GO TO ion 

INTEGFP FUNCTION MODLST t^X .PAPMStLQAy ,NUMO0) 
real WX (^*36^) *0AR’^S(4#366) 

gU«iSu??NF for EXTPFmElY ORV, CONglT 

Planting day ano« if needeo* recalculates 
CMECK^Vhf sponge for ORYNESS. return if it is MOIST. 

IF (!XX<6.L0AY) .GE. l.OOl RETURN 

DROUGHT Ml! 

1 » 0 
FLAG « 0 

NUMPO « NUMPn - I 
JO « LOAY - I 
100 JO « JO ♦ I 
CHECK PL AMT I MG OATF 
IF (PARM«;| I.JOI 
TF (PAPMSIl.JOl 


S ON THE LA 
LAST PLANT 




CRITERIA 
•LE. 2.0) 
.LF, 2.0) 


1 


NUMPD * I 


I ■ I 

NUmPO 

IF ENOUGH PRccTOTtATiON OCCUR*} HEFORE.A CRITICAL 
THFN WESFT THK DAY COUNTER AK-0 SET THE FLAG TO 
IF IFLAG ,e0. 0 .AMO. 

< (WXn.JO) .GE. 0.30 ,Awn, PARMS(4.J0) 
IWXn.JO) .GE. 0.7S .ANO. PARMSI4.J0) 

IF (FLAG .EO. 0 .AND. 

{ (RXO.JO) ,GF. 0.30 ,AN0. PARMS(4,JO) ,LT. 

‘ ~ 0.7S .AND. PARNS(4.J0) .LT. 

100 


A 

A 

1 

A 

A 


IF (I 


MOQLST « JO 
TTURN 


(WX(3.JO) .GE. 
.LE. 10) GO TO 


VALUE QF„6Q032R.„ 
THE fST RAINY OAY 

•LT. 600.00) .OR. 
.LT. 600.00))) I* 


600.00) ,OR. 

800.00) FLA6«J0 


REPRESS NT. A^T .1 V E„ O^A T E S 


RE1 
END 

»«*•«• OFTFpMlNE _ 

SUBPOUTINF ILASTC IFIRST.LASTfN.RX.P4RMS*LT0TALfNUHPD.M0L) 
REAL WX(6.36»^) .PA«MS(4»3A6) 

INTFGFR MnL(lOO) 

deteomimf last, the last planting date in the.string 
OETEPMInF M. the number of planting OATES IN THE STRING 
5 jn » IFIRST - I 
N ■ 0 

JO ■ JO ♦ 1 

IF (pARMS(l.jn) .LE. 2.0) N « N*l ^ 

IF(PARMSn .JD) .LE. 2.0) MOL (N^LTOT AL) «JO 
MN^LTOTaL) .EO. MUMPO) GO TO 20 
(PARMSn.JG) .LE. 2.0) GO TO 10 

- -- -- THE FIRST ( 

IFIRST ♦ I 


10 


If 

r.PST 
IF 

I . Jh»l) 


DAY IS WET. SO RESET 
(N .LT. 1) IFIRST * 
<N .LT. 1) GO TO S. 
WFT O/'Y— IGNORE |T 


DATE ANO START OVER 


OAYS— ENO 
I 


OF 


2.0) GO 
STRING. 


TO 10 
RESET 


THE JUUAN OATE ANO RETURN. 


MULTIPLE wet 
JO ■ JO •* 

?0 LAST * JD 
RETURN 

SUSPOUTINP SPl IT(IF.ME0.TL.N.nAYl.0AY2.0AY3) 

INTEGFP OAYl.OAY2.OAY3 ^ 

RETURms up to 0 OATFS TO REPOESFNT STRINGS OF PLANTING OATES 
values APF OAYl. OAY2. OAYl 
K,»S APE SUWSTPIMG LENGTHS 
IE*5 ARE FTBST OAYS IN SUBSTRINGS 
OAYl ■ 0 
DAY2 « 0 
OAY3 * 0 

IE (N ,GT. IS) GO TO 15 
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omtiluAi. ^ 
OF POOn QUAUTi 


ritn K0»v» A 


coMVFBSArioNAL ^ONim sysreM 


C STR 


c 

c 


i 


c 

c 

? 

8 

c 


!f« 


in 

,OMa -- f)N€ 

»M»i-n,ru 


8 


(M .r.T, <*,0 TO 

Mfi IS HP TO JO^OATS LO 
yi « IF . mosFS(N.|F. 
prTURH 

STPIMG IS ll-is OAYS LOMG 
TWO P£PPF«!FNTAT|''E 04VS 
10 K a N/> 

KK a N - 
IlF a IF * « 

OAVl a IF ♦ MnSF5(K.IF.MA0*(llF-lM 
0AV2 a IIF ♦ MOSeS<KK*tIFtMFO.IU 
RFTUPN 

STPIMG is 16* DAYS LOMO 
TMRFP RFPPFSFNjTATIVt PAYS 
15 K a N/3 
KK a N - Y 
KKMOn a Mftn(KK*2» 


PgPHESE^TATlVC DAY 


ONC sxTpa OAV — 
IF (KKMOO ,NF, 

. <><Kwon ,nf.. 

TWO FXTPA HAYS — 
TF IKKMOn .FO, 
F^tKKMOO ,E0, 
TF a IF * K 
IIP a IIF * 
JaYI a if * 

OAVS a 
0AY3 a 
RETURN 
END 


ASSIGN IT TO CENTER SUfjSTRINCi 
0) KK a K * I 
f>) KKK a K 
SHIFT TO LATFR OATES 
0) KK a K * 1 
0) KKK a K * I 


kk 

MOSESIK.IF.MED* (IIF- 

iF ♦ MOSES (KK* np*HEo. r 

IIF * MOSES (KK* 1 1 IF. MED 


•IlH* 


m 


.FUNCTION 


MosPs(N.n.MEn.iu 

OF DAYS, TO 


OETEOmINF number of days To median of string of length n 

EVEN STRINGS* MEDIAN IS SHIFTED TOWARD THE OVERALL MEDIAN PLANTING 




jPEMOP 


0.5 




FLOAT (Ni /a. - 
MO(HN.?> 

(IRFMOO.EO.O .ANn.^II.LE, 
MOSFS a MOSES ♦ T 


MED .AND. lL.LE.(MED«MOSESn 


»#•*•••« OUTPUT SUBROUTINES 

SUBPOUT INP ANSWFR(YR.ISEG.IPO.MPn,LPO.PLANT.PlNOEX*MDL*NUMPO> 
ALL CALLING oarAMETERS ARE INPUT PAPAMETEHS. THE 
SUBROUTINE PoODUCtS ONLY PPINTPf) OUTPUT— THE RESULTS OF THE MODEL 
YR^(YEAR). ISEG (SFGMENT NUMBER). ARE LOCATION IDENTIFIERS 
IPO. MPO. LPD. AMO ARRAYS PLANT ANO MOL ARE MODEL RESULTS 
PLANT contains THE PERRESENT AT I VF DAYS 
MOL contains all MOOELID PLANTING OATES . . 

PINOEX IS THP number OF DAYS WHICH REPRESENT THE PLANTING PERIOD 
(THP NUMPFR OF USPFUL DATES IN PLANT! 

NUMPD IS THE NUMBER OF MODELED PLANTING OATES, IN THE PLANTING 
PERIOD (TMF number OF OATES IN MOL FOR THIS SEGMENT! 


REAL PARMS(4.166! 

INTEGER MDL< lOO) *YR. PLANT (?R! 


WR^TEiT.Tnort! 




X! 

AN 

OF 


REPRE- 


^ _ YR. ISEG* IPD. MOD. LPD* (PLANT (n.laj, 

write year, segment number, initial planting date* me 

PLANTING date. LAST PLANTING DATE. AND THE RINOEX NUMBER 
C SENTATIVF OaTFS GFNERATED BV THE MODEL 
7000 F0PMAT(1?.1RI4! 

WOITF(7.7niO) (MDL(I) » ;*l .NUMPO! 

C WRITF the modeled PLANTING OATES (THERE APE NUMPO OF THEM* BUT 
C AT LPAST ??! 

7010 FORMAT (20T4! 

RETURN 

C 

SU8RDUTTNF TFST (WX .PARMS* ISEG! 

C TEST WRITES MftHTvuM aNO MINIMUM TEMPERATURES* PRECIPITATION. SPONGE. 
C SPONGE/PPFCIoITaTION variable* GD032R*4S« GOD* AND GD032R TO FILE «. 
ISEG IS AN Identifier, 

WX AND PARmA, AOg STaNDARO ARRAYS, 

PFAL WX (6.364) ,PARmS(4. 366! 

on To |a9o,isn 

WRITF (B.flOnO! ISEG.If (wX(J.I) *J*l*3! fWX(6«I! * (PARMS ( J* I ! * J*1 *4! 
10 CDNTIMUF 

8000 FDPMAT (2I4.RF7.2! 

RETURN 

END 
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r’ACE IS 
OF POOR QUAUTY 


PLDPvo foptoam ^ 


CONVKftSATIONAL MONITOR SYSTEM 


^ I L L 
SURPOOTfNP Pf 


PFAL 

tMTEGEO “ “ ' 
PFAL L AT 
PCAO HFAOFP r^Pr^ 


w X A P P A V •••••••••••• 

TPnanSFGt»*X»LAT.DIV*YP.ST*CPDI 


.VP*STfCPn*P(lA.> 


ATtfPD.OIV 
f l?Xf A4) 


20901 


AT 

THE 


' nf 


A," 


EAn the flag capo before oftoCEEOiNG TO 


9 

I? 


PFAn(?<i.prtnno.EMn«pp) iseg.st* 

20090 format ( f4.4X« r2f ;>X*FS.2f A2X«I 
00 10 »»'CAons«i.PS 
C PEAO PBFrTPrT^TION CAROS 

PFAO«Pn,9000NENO«4p) JOtPtIVR 
FoRMATIAX.n^ lX. lM4.4Xt rz) 

F (NCAPns .pfi. u yp«ivil 
F<JO.EO.-PP>GO TO 15 
0,5 ♦■1.16 
N0»J0.J-1 

IFIN0,5T,36A>r,0 TO 
THE EMO OF THF YEAP. 

TEMPFPATDRF OATA 
WX(1.MO»»P(J)*0.ftl 
5 CONTINliF 
COHTTNKE 
__ CONTI HHF 
C PEAO TEMPEPftTUPE OATA 
00 25 MCAPOSsIt^H 

oEA0t26.P000?.FN0*9Q> JO. ( (T(K.J> .K«l*2) .J"l.lO) 
F09HAT(6X.I3.PO!3» 

TF(JO.FO.-hA)PETUPN 
00 20 Jal.lO 
Nn«JO*J-l 

IF(N0.6T,36(«,>r,0 TO 25 
wx(i.Nn)«Tn.j)»j.o 
wx<2.Mn)«T(2»JJ*i.o 

cSSnAiT 

RFTUPN 
STOP 
ENO 

SUPPOMTINF SPIN (STN.YP.SPINIT. CAP) 

C GIVEN THF I 00 AT ION (STN) AMO YEAR (VP). READ FPOH 
C THF INITIAL VALHF FOR THE SPONGF. IfSUALLY. THIS 
J JpAST OAY OF THE PREVIOUS YEAR FOR THE LOCATION 


20002 


?0 

25 


99 


!l 




HE SPONGF IS IMTTIALIZED AT HAUF-CAPAC I T Y <CAP/2) 


THt 


19 


ON the 
NO VALUE. 


INIEGFO STN.YP.IO 
‘ 'THF OUTPUT 


parameter of THF INITIAL VALUE OF SPONGE. 


C SPINIT IS 

SPINIT ■ rAP/2,0 
C DEFAULT VALUF OF SPINIT IS HALF CAPACITY 
IF(YR.LT.O0)Pcf||tf.j 

C ONLY YEAPS after iPHO HAVE INITIAL VALUES. 

S PFA0(1P. 19000. EM0»10) 10. VAl UE 
19000 FORMAT nS.FS.?) 

IF no .EO. STN) SPINIT ■ VALUE 
C match CUOoF.NT 10 WITH STN TO GET PROPER VALUE OF SPINIT 
60 TO 5 
10 PFWlNO 19 
RETURN 
ENO 

SUBPOUT T NF SPONGE ( HX . C AP . SP IN I T ) 

C THE SPONGF moISTUoE VApTaBI.E. m.H. TRENCHARD. 

C INPUT STANOapo vx array OF metfoPOLOGICAL VARIABLES 
C SPONGF CAPACITY CAP. AA|0 STARTING VALUE OF SPONGE SPINIT. 

C OUTPUT OAlLY sponge VALUES IN WX(*.4) 

REAL WXf6.146) 
spongf*spinit 
no 5 jo*i 

C CALCULATE Daily VALUES 

TXawx ( I .JO) 

TM»WX{?.Ji) 
poMxwX (3. tO) 

TX-MAXriJM TFnpERATURE (F) .TN-minimum temperature IF). PPN-PRECIPITA* 
TION (IN). 

SOOf'OFsSPONGF.PPN-FP(TX.TN) *SP0NGF/<30.0*CAP) 

EP IS TPENCNAtiVS monthly PAN EVAPORATION FUNCT ION--OIV lOE BY 30 DAYS 
TO GET A OAiiY Value 

IF ( SPONGE. GT. CAP ) SOCMGEaCAP 
IF(SOONGF.LT.n.O) SPONGEsn.O 
5 WX (6.JD)sSP0NGE 
RETURN 
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ORIQIMAI- PAQ6 13 


pt novp rnpTPAM a 


CONVCPSATlnNAL MON I TOM SYSTEM 


C 

f 

i 


C 


fMO 

hlNCrnN r«(TY*TM> 

M.M.T»ENfiMAPn?s MONfHLy pam evapopatton ruMCfinN fP, - - 

INPUT MaxfMif*. nNO MiNlMtiM fPMMfPATlWPS <Ta ANO TmT In 0E6« F. 
Output pam EWAonaATloN in imch*.-s, wtu. r, 

USE, VAPOP function to calculate SATilPATlON VAPOW PPESSUkC 
f 275: ( Tx » -0 • p6^4#vapop < tn> 

^ IFIT»,LF.1P.0)Eo»0.0 

maximum TEMPF.MaTUME is rtELOM FPECZING, 

PIFTUMN 


pio 


«aTE». 

MM 


UMCTTOM *MrOO(T» 
^ NCHAPr ' - - 

INPUT 


M*H',f»CNCHAPOr< FUNCTION TO CALCULATE SATURATION VAPOR PPESSUPC OVER 
TFMP£P A TUPF T IN OEO. Ft OUTPUT VAPOR PRESSURE IN 


*AfTAinAi ' 


RPfURN 

END 


NASA-JSC 
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